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Abstract 


Let K be either the real, complex, or quaternion number system and let O(K) be 

the corresponding integers. Let X = ( Xl X n ) be a vector in K n The 

^torX has an integer relation if there exists a vector 171 = (iTli m ) S 

0(K) m * 0, such that m lXl + m 2 x 2 + . . . + m n x n = o’. In this paper we 
define the parameterized integer relation construction algorithm PSLQM, where 
the parameter Tcan be freely chosen in a certain interval. 

inning with an arbitrary vector X = (x-j , . . . , X n ) ^ K n , iterations of 
r oLQ(r) will produce lower bounds on the norm of any possible relation for X. 
Thus PSLQ(r) can be used to prove that there are no relations for X of norm less 
than a given size. Let M x be the smallest norm of any relation for X. For the real 

□ci C S P ! eX CaSC and CaCh fixed P arameter rin a certain interval, we prove that 
roLU(r) constructs a relation in less than 0(rT + n 2 log M x ) iterations 
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1. Introduction 


Let K be either the real, complex or quaternion number system and let 
O(K) be the corresponding system of integers (i.e., ordinary integers, Gauss- 
ian integers, or Hamiltonian integers, respectively). Let x = (xi, • • • ,-Cn) be 
a vector in K" . The vector x has an integer relation if there exists a vector 
m = ( mi ,. . . ,m„) € O(K)", m ^ 0, such that miXi + m 2 x 2 + ■ ■ . + m n x n = 0. 
In this paper we define the parameterized integer relation construction al- 
gorithm PSLQ(r). The parameter r can be freely chosen in the interval 
1 < r < p, where p is 2, v/2, or 1, depending on whether K is the real, 
complex or quaternion number system, respectively. We analyze PSLQ(r) 
for these three number systems. We describe in detail some efficient For- 
tran multiprecision computer implementations of PSLQ(r). We also present 
working Mathematical code for PSLQ(r) and, for comparison, some other 
relation finding algorithms from the literature. 

Beginning with an arbitrary vector x = (aq, • • • , x„) G K" , a finite number 
of iterations of PSLQ(r) will produce lower bounds on the norm of any 
possible relation for x. The computation of such a lower bound constitutes a 
proof that x has no integer relations whatsoever of norm less than this lower 
bound. Any finite computation can only prove that no small re lation exists . 

Let M x be the smallest norm of a relation for x. Let r = 1 / yl / P h V 7 > 
where p — 2 for the real number field and p = \/2 for the complex number 
field. For each fixed parameter r in the interval 1 < r < p, we prove 
in the real and complex case that PSLQ(r) constructs a relation in less 
than (") log r (7 n ~ 1 M x ) iterations. This shows that PSLQ(r) is “polynomial 
time” in the dimension and the number of bits of a smallest integer relation. 
Different r or 7 choices lead to different time and space requirements for the 

algorithm. 

For dimension n = 2 we prove that PSLQ(r) will construct a relation of 
smallest norm M x . We give examples in dimension n — 3 , for some r, for 
which PSLQ(r) does not construct a relation of smallest norm M x . However 
for any dimension n > 2, we do prove that any relation constructed by 
PSLQ(r) has norm less than or equal to 7" 2 M X . 

The “polynomial time” and “small norm proofs given here are straight- 
forward generalizations to the parameter r and to the complex numbers of 
the original “polynomial time” proofs that appear in Lagariaset al, [HJLS 89 ]. 
We show, however, that the algorithm of [HJLS 89 ] is distinct from any of 
these PSLQ(r) algorithms. 

PSLQ(r) was introduced by the authors [BaFe 91 ] in 1991 . PS refers to 
partial sums of squares, LQ to a lower trapezoidal orthogonal decomposition, 
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and (r) is a parameter defined as above. Since PSLQ(r) was introduced it 
has been used to discover numerous previously unknown identities among 
real numbers. One example is 



where L n {x) denotes the polylogarithm function Ylk xkfc ~ n - See [BaBG94] 
for details. Another example is the following formula for i r: 


00 l 


1 


;=o 


16 * 82 + 1 82 +4 Si + 5 82 + 6 


)• 


This remarkable series permits one to rapidly compute individual digits from 
the hexadecimal expansion of tv. See [BaBP95] for details. It was found 
by applying PSLQ(r) to the vector X = {X U X 2 , -- , AT 8 , 7 r) where Xj = 
Ylk>o l/(16 fc (8A: + j)). The smallest relation known, 


(4, 0,0, -2,-1, -1,0, 0,-1), 

yields the above “base 16” formula for n. A next smallest relation known, 

( 0 , 8 , 4 , 4 , 0 , 0 , - 1 , 0 , - 2 ), 

was subsequently discovered by Ferguson and this relation yields a similar 
base 16 formula for n. Together these two integral lattice relation vectors 
generate a two-dimensional lattice of relations of this “base 16” type. It is 
conjectured there are no further such relations outside this lattice. Note that 

(- 8 , 8 , 4 , 8 , 2 , 2 , - 1 , 0 , 0 ) 

is in this lattice, so evidently X 7 is integrally dependent upon Xi ,.Y 6 . 

Of course, a numerical discovery of a relation using PSLQ(r) does not 
constitute a rigorous proof of the relation. However, in the wake of this 
numerical evidence, proofs have subsequently been found for many of these 
relations, including the above formula for tv. See [BorBG95] and [BaBP95] 
for details. 

In the theoretical proofs in Section 2, 3, 4, and 5, we will assume exact 
arithmetic over the real numbers augmented by comparisons over the reals 
and the nearest integer function. 
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2. Lower Bounds on Integer Relations 

If K is the complex number field, then 2* denotes the complex conjugate 
of 2, i.e. if x = x + iy, then z* = x - iy. | • | denotes the complex absolute 
value, i.e. |r| 2 = 2*2 = 22* = £ 2 + y 2 • If A is a matrix or vector, then A * 
is the conjugate transpose of A. A unit in the complex number field is any 
element 2 such that \z\ = 1. For real 2, the conjugate operation is null, and 
2 is the usual absolute value. 

Similarly, if K is the quaternion number system, then 2* denotes the 
quaternion conjugate of 2, i.e. if 2 = x + yi+uj+vk , then 2 = x — yi—uj—vk. 
The quaternion absolute value or norm is similarly defined, so that \z\ = 

zz * _ z * z = X 2 + y 2 + u 2 + v 2 . Units and conjugates of matrices are defined 

analogously. 

If K is any of the above three number fields, two vectors x, y G K” are said 
to be orthogonal if xy* = 0. Let \A\ = (tr(A*A)) 1/2 denote the Frobenius 

norm of the matrix A, i.e., |A| = (£<>«<.>) 1/2 • An n x n matrix A is 
unitary if A* A = AA* = U{n, K) denotes the group of unitary matrices 
over K. An n x n matrix A is unimodular if det A is a unit. GL(n, O(K)) is 
the group of unimodular matrices with entries in the integers O(K). 

Definition Is (M r ). Assume x = G K" has norm |x| = 1. 

Define x x to be all vectors in K" orthogonal to x. Let O(K)" n x- 1 be the 
discrete lattice of integral relations for x. Define M x > 0 to be the smallest 
norm of any relation for x in this lattice. 

Definition 2: ( H x ). Assume x = ( x 1 , . . . , x „ ) G K n has norm |x| = 1. For 
1 < j < n define the partial sums 

s 2 = Xfcx£. 

j < n 

Given such a unit vector x define the n x (n - 1) lower trapezoidal matrix 
H x = (h,, } ) by 

{ 0 if 1 < i < J < n - 1 

Si+ifsi ifl<i=j<n-l 

-x*Xj/(sjS j+l ) ifl <j<i<n. 

Note that h itJ is scale invariant. 

Lemma 1. Let H x be the lower trapezoidal matrix defined above. Then 
(1) H*H X = I n - *.e., the columns of H x are orthogonal, 
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(n) \H X \ = v4n, 

(ill) XHj = 0. 


Proof. The columns can be proven orthogonal by considering the cases i = j 
and ^ < j separately. When i = j the inner product is 




ti . V _ s *+i , x,x* ^ 

2 ^ shL, + 


i<k<n *i *i*U 

_ f? - flf-l ^ XiX* 

When i < j the inner product is 


+ 

•S" s* 


= 1. 


Sy^jXjXj xJXjXj^ fjc 

SjS,S, + 1 ' ^ S t S l+1 SjS J + l 


_ _ fj+ 1^,-gj xjxj_ y- * _ 

•S,Si +1 s J s J+1 E ~ 

1 1 J<k<n 

Item (i) shows that HfH x = / n _j which has trace n - 1 so \H X \ = ,/n - 1 
To prove (iii), fix 1 < j < n - 1, then ' 

E ****,,• = fi£ ^ ±i - y wfk x i' = 

TTT. S , ^ a e ■ . . 


1 <k<r 


„ S }Sj + 1 


J<k<n 


x j s j + 1 


s j s j + 1 


= 0. 


□ 


p F ° r a um * uec<or ^K n de/me P x = tf x // x *. 7%en P x satires; 

17 — *X i 

(n) P x = I n -x*x , 

(m) Pi = P x , 

M \Pr\ = 

77 -Pr^* = 2 * /or any 2 £ x -1 - 

( 01 / Pj.ni* = m* /or any relation m 6 Q(K)" for x. 

Proof Item (i) follows from P X P X * = (tf x tf;)\ To prove (ii) note that 
from Lemma 1 (an), H x is an n x (n - 1) rank n - 1 matrix whose columns 
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transposed form an orthononnal basis for I 1 . Defining U - lH r \.r’h an 
„ x n unitary matrix, we have UV = H,H‘ + r*x = /„. To prove (...) note 

that „ _ „ 

p 2 x = (/„ - x*x) 2 = II - 2/ n x*x + X (xx )x - Px. 

To prove (iv) note that |P,| 2 = tr(P,*P,) = trP x = tr H' t H, = n - L Item 
(vi) follows from (v) which follows from (n) and the associativity (x x)x. 

X * ( X 2 * ) . □ 

Theorem 1. Let x ± 0 € K n . Suppose that for any relation m of x and 
for any matrix A € GL(n, O(K)) there exists a unitary matrix Q 6 U(n - 1) 
such that H = AH X Q is lower trapezoidal and all of the diagonal elements 

of H, hjj ^ 0. Then 

= min tt — r < l m l- 

maxi<_,<n-i \hj,j\ \ h jj\ 


Proof Let m be any relation for x. By the hypothesis, there exists a uni- 
tary matrix Q 6 U(n - 1) such that H = AH X Q is lower trapezoidal (this 
is equivalent to QR factorization). There is an n x n - 1 matrix c T with 
diagonal ones and an n - 1 x n - 1 diagonal matrix D where H - TD with 
diagonal entries h jtj ± 0, 1 < j < n - 1 from the hypothesis. On the other 
hand, AP X = HQ*H *, from the definition of P x in Lemma 2. The equation 
_ TDQ*H* x gives a decomposition of AP X into the product of a lower 
trapezoidal matrix T with diagonal l’s, an invertible diagonal matrix D with 
diagonal h' s, and an n - 1 x n matrix Q*H* X with orthonormal rovvs since 
Q+H * H x Q = Q*In-\Q = /n-1 by Lemma 1. So the norm of the j - th row ot 

D %voi Lemma 2, part (vi), m* = P,m\ so that Am* = AP t m\ From 
the above decomposition of AP X = TDQ*H* X , we have Am = AP x m - 
TD{Q*H*)m*. Let Qhj be the j-th row of Q* H* and let Aj be the j - th 

row of A. Then 

A } m * = hj,jQ H ,jrn* +Y, t ^ khk ' kQH ' kTn *' 


Since A is invertible, Am * f 0. Let j be the least j for which A } m* ± 0 so 
that A k m* = 0 for k < j. Then the k < j rows of TDQ* HJm*axe zero, and 
since T is lower trapezoidal by recursion, the A:-th rows of Q H x m are a so 
zero. With this least choice of j then A } m* = h h) Q } m*. Therefore, from 

A € GL{n, 0(K)), 

1 < |A ; m*| < I hjjQ H,j m *\ — I’ 
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because Qh,j is a unit vector. □ 

Comment on Theorem 1. Theorem 1 suggests a strategy to construct 
a relation finding algorithm: Find a way to reduce the norm of the matrix 
H x by multiplication by some unimodular .4 on the left. The inequality 
of Theorem 1 offers an increasing lower bound on the size of any possible 
relation. Theorem 1 can be used with any algorithm that produces any 

GI(n, 0(K)) matrices. Any GL(n,Q(K)) matrix .4 whatsoever can be put 
into Theorem 1. 

Definition3: (Hermite reduction). Let H be a lower trapezoidal matrix, 
with h tJ = 0 if j > i and h h} ^ 0. Define the matrix D = (d, } ) £ 

GL(n,Q(K)) recursively as follows. For fixed i, decrement j from n to 1, 
setting 

( 0 if* < J 

d ttJ = s 1 if i — j 

[ nint((- J2j<k<i d i.khk,j)/h jtj ) if j < *, 

We will say that DH is the Hermite reduction of H and we will say that 
D is the reducing matrix of H. The function nint denotes a nearest integer 
function, e.g., nint(t) = [f + 1/2J . This definition of nint can be extended to 
each coordinate for complex or quaternion arguments. 

Definition 4: (Modified Hermite reduction). With the same notation 
as in Definition 3, set £> = /„. For i from 2 to n, and for j from i - 1 to 
1 (step -1), set q = mn then for k from 1 to j replace h l[k by 
ht.k - qhj,k, and for k from 1 to n replace d i k by d hk - qd J>k . 

Lemma 3. For a lower triangular matrix H with h,j = 0 if j > i and 
h } ,) 0, Hermite reduction is equivalent to modified Hermite reduction. 

Comment. This variation can be found in [BergSO] and later in [LLjL82], 
This recursion replaces the input H with DH while developing the left mul- 
tiplying reduction matrix D. 

Lemma 4. There exists a constant p K = p > 1, with the property that 
the entries of the Hermite reduced matrix H' — (h\ ) = DH satisfy the 
inequality “' J ' 

Wul ^ \ h U/p = \h t Alp 

for all k>i. The constant p = 2 for the real case, p = y/2 for the complex 
case, and p = 1 for the quaternion case. 

Proof. This follows from the definitions of the nint function, Hermite reduc- 
tion, and the fact that |z - nint(z)| < v'dim» K/2 for z e K. □ 
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3. Statement of the Algorithm PSLQ(r) 

Definition 5: (The parameters 7 and r). Fix the real number 7 > 2 /VZ 
or 7 > \/2 or 7 = 00 for the real, complex, and quaternion cases respectively. 
In terms of this 7, define the real number r by 

1/r 2 = 1/p 2 + 1/7 2 , 

where p is defined as in Lemma 4. For the proof of Theorem 2, we will 
require that 1 < r and that r < p; clearly these conditions are satisfied in 
the real and complex cases. In the quaternion case r = 1 and p = 1. 

For the proofs that follow assume K is real or complex, not quaternion. 
Note however that the statement of the algorithm is valid for the quaternions. 

Initial conditions: Given the input unit vector x G K n , set H = H x where 
H x is defined as above. Set the n x n matrices A and B to the identity 
Perform Hermite reduction on H, producing D G GL{n,Q{ K)). Replace x 
by 1P -1 , H by DH , A by DA, B by BD 

One four-step iteration: 

Step 1: Exchange 

Let H = where h l<} is the i-th row, j-th column entry of H. Let 

Q = h r< r , (3 = h r+ i.r, A = fi r+ l,r+l, $ = \[W* + AAV 

Choose an integer r such that 7 r |^r,r| > f° r all 1 < 1 < n - 1. Define 

the permutation matrix R to be the identity matrix with the r and r + 1 
rows exchanged. Replace x by xR , H by RH, A by RA, and B by BR. 

Step 2: Corner 

At this point the updated matrix H may not be lower trapezoidal since 
A may not be zero. If r < n - 1 replace H by HQ where Q is the unitary 
n-lxn-1 matrix Q = (q tJ ) € U{n - 1, K) defined by 

f3* 1 8 if 1 — r,j — r 

— A /S if i = r,j = r + 1 

A */8 if i = r + l,j = r 

3jS if i = r + l,j = r + 1 

1 if i = j ± r or i = j ± r + 1 

0 otherwise. 
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where the a,/?, A, are defined in Step 1. If r = n - 1 then H is unchanged. 

Step 3: Reduction 

Perform Hermite reduction on H, producing D € GL(n, 0(K)). Replace 
x by xD 1 , H by DH. A by DA, B by BD~' . " P 

Step 4: Termination 

Terminate the algorithm if Xj = 0 for some 1 < j < n or if h, , = 0 for 
some 1 < i < n - 1. ’* 


4. Number of Iterations of PSLQ(r) 

of PSLO ( |lr H ' a’ a /M S 7,r be the result after exac “y k iterations 
of PSLQ. Let a = KAk) and /J = h r+I _ r (k). These definitions of q and 3 

are consonant with those of Step 2. Because H is Hermite reduced in Step 

< “l/e For C < n - 1 set A = h r+Ur+l{k ) and define 
t by t - y/00 + AA /for |. From this definition of t we have 

I'M < M t. 

From the Step 1 Exchange, 0 < |A| < |a|/ 7 . It follows that 


t = y/W* + AA*/|q| < y/l/p* + l/ 7 2 _ 




as in Definition 5 For this proof we will require that t < 1 < r , clearly 
satisfied in the real and complex cases. 

Lemma 5. If h hJ {k) = 0 for some 1 < j < n - 1 and no smaller k, then 
3 l and a relation for x must appear as a column of the matrix B. 

Proof. (Alyson Reeves) First we show that h h} = 0 implies that j = n- 1 
Consider the matrix H(k - 1), the end result of the k - 1-th iteration. By 
he hypothesis on k we know that no diagonal elements in H(k - 1) are zero 

kn pwT’ th ^ r ^ bout to be chosen Step 1 of the *-th iteration, we 
now that hr, r (k - 1) ^ 0 and that h r+1>r+1 (k - 1) ± 0. Now, suppose the 
r chosen in Step 1 is not n - 1. Let 

a 0 

0 A 

rln dVfr Aft X l ‘ ’ , C ° nsistin « of the ' “d r + 1 rows of columns 

+ 1. After Step 1 has been performed this submatrix becomes 


:) 
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At Step 2, we post-multiply the matrix by the unitary sub-matrix of Q 

(F/s -VA 

\X*/8 0 /*)' 


where 5 = sjW + A A*. The result is the matrix 


( * 0 ) 

\ap*/6 -aX/SJ’ 


Since A and a are not zero (they were diagonal elements of H(k - 1)), we 
know that and -aX/8, the two diagonal elements in the matrix, are also 
not zero. Note that since the rest of Q is the identity matrix none of the 
other diagonal elements is affected by the multiplication. Thus, at the en o 
Step 2, all diagonal elements are non-zero. Since Hermite reduction doesn t 
introduce any new zeros on the diagonal, the end result of the fc-th iteration 
has all non-zero diagonal elements. But this contradicts the hypothesis on 
k and our assumption that r < n - 1 was false. Note that for r = n - 1 
in order to have h„- i,»_i(fc) = 0, we must have h n ,„-i(k - 1) - 0 and 

h n -i,n-l{k — 1 ) 7 ^ 0 - f , • 

Next we show that a relation for x must appear as a column of the matri 
B. By Lemma 1, xH x = 0. BA = /„ implies 0 = xBAH x = xBAH x Q = 
xBH(k - 1), where Q is an appropriate unitary n - 1 x n - 1 matrix. Let 

2 = xB. The above gives 


(0,... ,0) = xBH(k - 1) = zH{k - 1) = (... ,Zn-lh n -l,n-l(k - !))• 

Since - 1) ? 0 then z n -i = 0. Hence the n - 1-th column of B 

is a relation for x. □ 

Lemma 6. At any k-th iteration of the algorithm the diagonal entries of 
H(k) satisfy the inequality |h ti ,(A:)| < 1. 

Proof. We follow the a,/?, A definitions of the proof of Lemma 5 and use 
induction. For k = 1 the diagonal entries of H{k ) are those of H x and 
s +i < sj < 1 gives the required inequality. Assume that the inequality also 
holds"up J to k - 1. The diagonal entries of H(k) are equal to those of H(k - 1) 
except for row r where Step 1 Exchange occurs. When r = n - 1, after the 
exchange, the r-th diagonal element is (3. But \,3\ < \<*\/p< 1 because p > 1 
and | Q | < 1 by induction. When r < n - 1, after the exchange the r-th 
diagonal element is <5. But \S\ = \a\t < 1 since f < 1 and |a| < 1- The 
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r + 1_th dia gonal element of H is -a\/S (as in the proof of Lemma 5) so 
that | - a\/ 8 \ = \\\jt < |o| because |A| 2 < |A| 2 + |ii | 2 and |A| < |a|t. □ 

We show that every iteration of PSLQ causes a geometric monotonic in- 
crease in a certain function Il(k) which is roughly the product of all the 
principal minors of the matrix H(k). If a relation for x exists, this prod- 
uct will be bounded above and below. Assume x has some relation and as 
usual let M x denote the norm of a smallest relation for x. We will need the 
following technical lemma in the proof of Lemma 9. 

Lemma 7 . Consider the quotient 

q(A n<) _ min{g,t}.min{.4.1} 

min{£?, 1} • min{A, <} 

Suppose that the four positive real numbers A,B, l.t satisfy the three inequal- 
ities 


Then, 


A> B, A>t , 1 > t. 

q{A,B,t) > 1. 


Proof. Of the 16 possible choices in the min’s, the inequality A > t removes 
8, A > B removes 2, and 1 > t removes 1 leaving 5. These five are 
A > B > 1 > t with quotient tj 1 • 1 jt = 1, 

A > 1 > B > t with quotient t/B ■ l/t > t/l ■ l/t = 1, 

1 > A > B > t with quotient t/B ■ Ajt = A/ B > 1, 

1 > A > t > B with quotient B/B ■ A/t = Ajt > 1. 

^ > 1 > t > B with quotient B/B • l/t = l/t > 1. □ 

Lemma 8. For a, 7, M x as above , 


l n ~ 2 M x \a\ > 1. 

Proof. By the choice of r in Step 1 Exchange, we have 7 r |o| > y'l h I for 
an y Ji 1 < J < n — 1 , which implies 

7 l\hj,j\ > 7 r /| hjj\ > 7 J /|or( > 7 1 /|a|, 

for all j including that j 0 for which M x > lf\h ioijo \ from Theorem 1. Thus 
7 2 M X > 1/|q| and 7'- 2 M,|a| > 1 □ 


Definition 6: (The II function). Recall r = ^/Tj^TVPf. Define 
n(A:)= [J min{ 7 n - 1 M J ,l/|h^(A;)|}"->. 


li 


Lemma 9. 

(i) 


For any k > 1 we have 

(7 b- 1 m,)(>) > n(fc) > 1, 


(ii) 

n (fc) > rn(fc-i). 

Proof. For the fc’ s so far, hjj(k) ± 0 for all 1 < j < n - 1. M x > l and 
1/| hjj{k)\ > 1 by Lemma 6. This gives 

min{Mj;, l/|hj,j(fc)l} > T 

for all 1 < j < n — 1, which implies the right hand inequality of (i). On 
the other hand, it is always the case that M x > min{M x , l/|h>,;(fc)|}) which 
together with the fact that (") = n - 1 + • ■ • + 2 + 1 and that 7 > 1 gives 
the left hand inequality of (i). 

The proof of part (ii) is more involved. Let r be given by the btep 1 
Exchange of PSLQ. Recall the definitions of the tw o successiv e diagonal 
elements a, A, the single off diagonal element (3, t = \/0(3* + AA*"/|a| in the 
Step 2 (Corner development) of the unitary matrix in terms of / 3 and A. 

Suppose that r < n - 1. Then only two diagonal elements change. These 
correspond to the 2x2 submatrix of H 

(;:) 

which after a single iteration becomes 

(8 0 \ 

\a/3*/5 -a\/8 J 

But |<S| = \a\t so that the absolute values of the of the a, A diagonal elements 
are replaced by the absolute values of the 8,—aX/8 diagonal elements. All 
the factors of II(fc) are the same except these two so that 

n(fc) ( m inb n - 1 iV/x,l/(M*)} y~ r ( min H-;^V|Ali y 
Il(fc — 1) \ min{ 7 n_1 M x , l/|a|} / \min{ 7 " 1 M x , 1 / 1 A | } / 

Set 

A = 7 n-1 M r |a|f and B — ~< n Af x |A|, 
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so that 


— n(^~) = ( /min {B,t} minimi, 1 } 

n(* — 1 ) \ min{T, t} 



) n — r— 1 

' 

We now show that the assumptions for Lemma 7 hold. Note that 1 > t 
by the definition of <; also, .4 > B since |a|t > |A|. By Lemma 8 we have 
A > *7 > By Lemma 7 we have 

H( fc ) > min {y4, 1} ^ l 


> - > r. 
f ~ 


n(k - 1 ) min{/l, t} 

Now suppose that r = „ - 1 . By Step 3 Reduction, under one iteration 
the absolute value of the last diagonal element a is less than |o|p. All the 
factors of II(A-) except the last are the same so that 


jjjjO < min{ 7 n l M x ,\/{\a\p)} 


H(^ - 1) min{ 7 n_1 A/ x , l/|o|} 


But we always have 7 n ~ 2 M x 


mm{A,t / p) 
min{/l,<} 


M > 1, so if A > t/p > t 

H(Ar) 


n(fc- 1) 

By Lemma 8,A>t~/>t. If t < A < t/p then 

Il(Jfc) 

n(4- - 1) - A ! t - 7 - r - 

Thus for r < n - 1, II(Jb) > rll(fc - 1). □ 

Theorem 2. Assume real or complex numbers, n > 2 r > 1 and that 

°ff X %f n h Zf^\ int TI rdatl ° nS • Let Mx be the leaH n0rm °f relations 
for x. Then PSLQ(t) will find some integer relation for x in no more than 

'n\ lQ g (7 n ~ l M x ) 

2j log r 

iterations . 

-"I? /nfuT "I have ,f° ne k iterati ° ns . ‘hen from Lemma 6 and Lemma 
*’ f 0 and not ' Jl 17,1(01 < 1/Mr for l < k. By Lemma 6, 11(0) > 1 

and by Lemma 7, IT(A-) > T k so that 

(7 , *“ 1 iW x )^s) > r k 

Taking natural logarithms of both sides of this inequality gives 


Q)l°g(7" l M x ) > klog 


□ 
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Corollary 2. Let K be the real numbers R or the complex numbers C. Fix 
n > 1 and assume given a unit n-tuple x € K” which has a relation »,£ 
OIK)” of least norm M, . Then there exist, a 7 such that the algorithm 
PSLQ(t) will construct some 0(K)“ relation for X m no more than 

2 • (dimiR K) ■ (n 3 + n 2 log A/*) 


iterations. 

Proof. Let 7 = 2. Then for either K, r > 1, specifically, 1/ log r < 4dimR K. 
PSLQ(r) takes O(n) exact arithmetic operations per iteration, so in is 
sense finds relations in “polynomial time” 0(n + n logM r ). 

5. Upper Bounds on Integer Relations 

We compare the relation found by PSLQ to a shortest possible relation. 

Lemma 9. Suppose m ,s the relation found on the k + Let iteration so that 
i(fc 4- 1) = hnm-i(h) = 0 ani hn-l.n- l(*0 + 0- Then 

|m| = l/|/ln-l,n-l(^)l- 

Proof At this iteration we have developed the matrix A 6 GL(n, 0(K)) 
where the (n - l)-st column of A -1 by Lemma 5 is m and the vector Am - 
t hi i its only non-zero entry a 1 in the (n - l)-.t positron. Srnce 
AP = TDO Qm * = D-'T'Am*, where T f is the generalized inverse of T 
and D is a diagonal matrix with last entry /»„-!, which is also the 
last entry of D-'T'Am*. Because Q is unitary \Qm*\ = \m |. U 
Theorem 3. Let Mr be the smallest possible norm of any relation for x. 
Let m be any relation found by PSLQ(r). For all 7 > yffi M ™ ctor * 
and for all 7 > y/2 for complex vectors 

|m| < 7 n ~ 2 M x . 

Proof. Assume we are at the fc-th step of PSLQ where a Step 1 Exchange 
r = n _ i W as made with 7^ 0 and h n - i,n-i( + ) 

7 n_1 |fin-l,n-l(^)l — l J \ h }A k )\ 

for all 1 < j < n - 2 by the choice of r. Hence, by Theorem 1 and Lemma 8 

Mr > 1/ max |h,,i(fc}| > l 2 ~ n l\h n -i.n-Ak)\ = 7 2 ” n |^l- D 
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Comment on Theorem 3. For n = 2 , Theorem 3 proves that any relation 

c j m ^ ) founcl has norm M = iV/ x- In other words, PSLQ(r) 

hnds a shortest relation. For real numbers this corresponds to the case of 
the Euclidean algorithm, [Euclid, Book X], [Fowl 79 ], [Knut81], For complex 
numbers this corresponds to the case of an algorithm in [Schm 75 ], 

For rc = 3, let x = (113,343,311). This vector has a shortest relation 
— (/,— 15, 14) with the shortest norm |m x | = M x = 21.6794 This 
can be verified directly, cf., [Kann 88 ], [PoZa84], [Cohe 93 ]. ' On the other 

hand for r = 1 . 0000 ..., 7 = 1.1547..., PSLQ(r) in iteration 6 produces 
the relation m, = (24, - 7 ,- 1 ). Indeed 


M x < |mi| = 25.0199 . . . < 7 M x = 25.0333 .... 

This relation appears from a zero in the second coordinate of the rT 1 
vector. Continuing to iteration 8 gives the relations appearing from the first 
and second coordinates of the current xA s 1 vector, m 2 = (-17 -8 15) and 
m 3 = (41,1,-16) of norms 24.0416... and 44.0227..., respeitivdy The 
vector m 2 has smaller but not smallest norm. Continuing to iterations 9 
a j_i 10 f 1VCS the relatlons appearing from the first and second coordinates of 
X ™ 9 of ( 7 ^ -15, 14) and m 2 = (-17, -8, 15), so a shortest vector m 4 

was eventually found. In iteration 11 the h 2 , 2 (ll) = 0 condition appeared for 
the first time giving the relation m 5 = (- 10 , -23, 29) of norm 38.3405 .... 

This example is instructive in that various choices of the parameter r give 
different outputs. The “legal” r are such that 1 < r < 2 , although the 
P SL Q(r) sometimes works for “illegal” r outside of this interval. For the 

Tin T ' ,T L1, lteration 6 y ields m u S yields m 2 ,m 3 , 9 yields m 4 ,m 2 , 
and 10 yields m 5 . On the other hand, for r = 1 . 8 , iterations 4 , 5,6 all yield 
only the shortest length relation m 4 . For the “illegal” r below 0.7 and above 
2 T the algorithm cycles indefinitely. The end point r = 1.0 gives essentially 
the same outputs as r = 1 . 1 . The other end point r = 2.0 yields two new 
relations m 6 - ( 1 , -91, 100) and m 7 = (0, -311,343) of norms 135.2109 
and 463.0010. . . , respectively. 


6. Multiple Relations. 

A given unit vector -r € K n may have 0, 1, 2, or up to n - 1 relations. 
Unce a relation has been constructed, one of the coordinates of xB for the 
appropriate B e GL(n, Q(K)) will be zero, and the corresponding column of 
B will be a relation. The remaining n - 1 coordinates can be used to form 
a new unit vector in ye K". Apply PSLQ(r) to this </. Any second relation 
so found will be integrally independent from the first and can be referred 
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back to the original x. In this way as many as n - 1 integrally independen 
relations for x can be constructed. We omit here the tangent discussion of 
using classical lattice reduction techniques to find integer relations; this is t e 
case for the Recognize[ ] function in Mathematical which calls the function 
LatticeReducef ], cf„ [Cohe93], [CJOS83], [LaLS82], Lattice reduction there 
applies typically only to integer relations for integer vectors. Integer relation 
finding here is directed specifically at integer or Gaussian integer relations 
for real or complex number vectors. 


7. Variations of PSLQ(r). 

The algorithm PSLQ(r) as stated may be performed for various “illegal” 
r or “illegal” 7 , and under these circumstances will find relations for some x 
vectors. This can happen for 7 < * the real case, for 7 < ^2 in the 

complex case, and for 7 < °° »n the quaternion case, so that r < 1 and the 
conclusions of Theorem 2 or Theorem 3 make no sense or have no apparent 
content. The reason for this apparent anomaly is that for a specific n-tuple x 
the actual field or division ring constant p bound in Lemma 4 is not universal 
and could depend upon an input vector x. Say Pl gives a bound such as that 
of Lemma 4 for some special x or collection of them. Then there may e an 

“illegal” 7 so that r, = 1 /vW + W > L For such x one CO " ld eX P e , Ct 
to see some relation emerge before the number of iterations indicated by 

Theorem 2 for this r x = r . a , 

On the other hand, it is possible to use the real PSLQ(r) algorithm to find 
complex and quaternion relations at the expense of doubling and quadrupling 
the dimension. For example, suppose z = x + yi + uj + vk is a vector in H 
with vector components € R* Suppose the corresponding relation 

j s m = a + bi + cj + dk which is a lattice point in W" with integral vector 
components a,b,c,d € Z» Then 2 m* = 0 implies four integer relations 
among the interlaced and suitably sign changed coordinates of 2 . For the 
first set Zi^nK-r, - Mi - W " Mi) = 0 and ° nC C& ^ * 

PSLQ(r) to the real 4n-tuple (. . . , x 7 , y } , u v Vj, . . . )■ There are ree ^o ers 
which are similar. A relation for * will be in the intersection of the four 
associated lattices. Alternatively, one can give a PSLQ(r) algorithm along 
the lines of [HJLS89, Section 5. Finding simultaneous integer relations J. 


8 . Computer Implementation of PSLQ(r) 

The PSLQ(r) algorithm can be implemented using ordinary floating point 
arithmetic on a computer. Using double precision (be., 64-bit) arithmetic, 
relations of two or three digits in size can be recovered for n up to five or 
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so. Beyond this level, precision is quickly exhausted, and recovered relations 
and norm bounds are meaningless. Thus a serious implementation of PSLQ 
(or any other integer relation algorithm for real numbers) must employ some 
form of multiprecision arithmetic. The authors employed the MPFUN mul- 
tiprecision translator and computation package. The Fortran-77 version of 
this software is described in [Bail93], and the newer Fortran-90 version is 
described in [Bail95], A C++ translator that employs these routines is also 
now available. Alternatively, one may employ the multiprecision facilities of 
symbolic math software packages, such as Maple, Pari or Mathematica T,v/ . 

The descriptions presented here of computer implementation of PSLQ(r) 
are for the case of the real number system. Extensions to the case of the 
complex and quaternions number systems are straightforward, provided one’s 
multiprecision system supports these datatypes. 

One key to an efficient implementation is to utilize a simplified version of 
Hermite reduction and the associated update. As noted in Lemma 3 above, 
Hermite reduction can be done more efficiently by a triply nested loop. In 
fact, the update operations associated with Hermite reduction (updating 
x, H, A and B) can also be done in a loop of this form. Further, if these 
updates are done in this manner, then it is not necessary to compute the D 
matrix. This simplified scheme is as follows. In the initialization step, Her- 
mite reduction and the subsequent updates are replaced with the following: 

For i from 2 to n, for j from i - 1 to 1 (step - 1 ), set t = nint {h ii} /h } ,) and 
replace x } by xj - tx t - then for k from 1 to j replace h hk by h,, k - th Jjk \ for 
k from 1 to n replace a i k by a hk - ta jk and replace b k] by b(k,j) + tb(k, i). 

Step 3 is also replaced with this, except i is incremented from r + lton and 
j is decremented from min {* - l,r + 1} to 1. Here r denotes the row index 
selected in Step 1 . These more restrictive limits on i and j merely reflect the 
fact that t = 0 outside these limits. 

Obviously in a computer implementation some care must be taken in 
testing for zero. This is typically done by checking that the absolute value of 
the tested value is less than the “epsilon” appropriate for the level of numeric 
precision being used. Also, a run should be terminated if any entry of the 
A matrix exceeds the level of numeric precision being used (so that these 
integer values can no longer be represented exactly). 

The level of working precision required for PSLQ is generally only a few 
digits greater than the accuracy of the input x vector. Along this line, if 
one wants to recover (or to exclude) relations of size d digits, then the input 
data must be specified to at least nd digits in order to obtain numerically 
meaningful results. The significance of a recovered result can be measured by 
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noting the ratio between the multiprecision epsilon and the largest entry of 
the updated x vector when a relation is recovered. If this ratio is very small, 
such as lO -40 , then one can be fairly certain that the relation produced y 
PSLQ is a real relation. But if this ratio is only a few orders of magnitude 
below unity, then the result is suspect, and higher accuracy in the input 
data, as well as correspondingly higher working precision, is required. 

The above implementation is satisfactory for most applications. For more 
demanding applications, a “two-level” implementation is significantly faster. 
In a two-level implementation, most operations are performed in or mary 
double precision arithmetic, with occasional updates of multiprecision .arrays 
using multiprecision arithmetic. This two-level scheme can be described as 
follows. Here the prime notation is used to denote double precision approx- 
imations to multiple precision values. 

To initialize, perform the initialization step as described above using lull 
precision. Then perform a “double precision initialization . (1) set x ~~ 
xl max, , k.l and set H' = H; (2) perform a LQ decomposition on H , using 
double precision arithmetic, setting H' to be the lower triangular part; (3) 


set A' — B' = In- / ill A' A 

PSLQ iterations are then performed as above on the arrays x ,M , A an 

B' using double precision arithmetic. Some care must be taken to insure 
numerical accuracy in these iterations. Obviously these iterations must be 
halted before entries in A' grow so large (9 x 10 on IEEE systems) a 
they cannot be exactly represented as double precision values. In the authors 
implementation, double precision iterations are halted when the largest entry 
of A' exceeds 10 10 . Tests for zero in these iterations must reflect the accuracy 
of double precision arithmetic — the authors used an “epsilon” of 10 ^ 

here. As an additional measure to insure numerical integrity, the authors 
code aborts the modified Hermite reduction procedure (and restores arrays 
to their previous values) if the multiplier q exceeds 10 . 

When the double precision iterations are halted, either due to large entries 
in A', or to a tentative zero in x' or H\ it is necessary to perform a “mul- 
tiprecision update”: (1) replace A by A' A, replace B by BB , replace H by 
A'H, and replace x by xB (2) check for zero entries in x, using the mul i 
precision epsilon. If no zeroes are found, then a double precision initialization 
is performed, followed by more double precision PSLQ iterations. 

One detail has been omitted here. In some cases, the entries of^the up- 
dated x vector have such a large dynamic range (greater than 10 m the 
authors’ implementation) that when converted to double precision, additions 
and subtractions would produce results of questionable reliability. In these 
cases it is necessary to perform PSLQ iterations on the mult, precision ar- 
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rays, using multiprecision arithmetic, for a number of iterations until this 
large dynamic range is eliminated. If this situation is encountered on any 
iteration other than the very first, a multiprecision LQ decomposition of H 
must be performed prior to performing these multiprecision iterations (so 
that the H array contains the same entries as the H array defined in the 
PSLQ algorithm statement). 

9. Summary of the Literature 

The problem of finding integer relations among sets of rational and real 
numbers is quite old. When n = 2 this problem can be solved for ratio- 
nal by the first Euclidean algorithm in Euclid, Book VII, and for reals by 
the second Euclidean algorithm given in Book X, cf., [Ivnut81], [Cohe93], 
[Shim94j. Generalizations of this algorithm to higher real dimensions were 
proposed without proof by many authors, including Jacobi [Jal868], Her- 
mite [Hel850], Poincare [Pol882], Perron [Perr07], Brun [Brunl9, Brun57l 
and Szekeres [Szek70]. 

The first integer relation finding algorithm with proofs for the case of 
real numbers was discovered in 1977 by Ferguson and Forcade, [FeFo79, 
FeFo82], These algorithms were shown to be polynomial time in the loga- 
rithm of the size of a smallest relation. They were not shown to be polyno- 
mial in the dimension. Since then, other algorithms for finding relations for 
real vectors have appeared in [Berg80], [Ferg86]. [Bail88] reports on a com- 
puter implementation of [Ferg86]. The sequence including [HJLS89] (HJLS), 
[BaFe91] and [ArFe93] (PSLQ), [BaBG94] (a concise statement of PSLQ), 
and [RoSc95] (a stable variation of HJLS) will be discussed below. 

These algorithms all depend upon an orthogonal decomposition of some 
sort. See [GoVL90], for a list of various orthogonalization algorithms and 
their differences. PSLQ is of the QR type. HJLS follows the lattice reduction 
work of [LLjL82], which is classical Gram-Schmidt type, cf. [Cohe93], This 
difference may explain some of the differences observed between PSLO and 
HJLS, cf. [BaFe91]. H 

Rigorous proofs that the algorithm under investigation must find a relation 
if one exists appeared in [FeFo79, Berg80, FeFo82, Ferg86], All of these 
proofs gave a linear bound in the logarithm of the size of a relation, but 
were not known to be polynomial in the dimension. [Berg80] and [Ferg86] 
had unsatisfactory proofs in the sense that they were shown to be at worst 
exponential in the dimension rather than polynomial in the dimension. This 
unsatisfactory state of affairs was resolved affirmatively with the proofs that 
appealed in [HJLS89] for the “small integer relation algorithm”. We will 
refer to this “small integer relation algorithm” as HJLS, as stated in [HJLS89. 
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Section 3] as a reflection of that in [Berg80, Section 3]. In fact, this proof 
in [HJLS89] was the first appearance in the literature of a polynomial time 
bound for a relation finding algorithm, polynomial in both dimension and 

logarithm of relation size. ... . , 

This important progress was made when [HJLS89] combined two in e- 
pendent streams of research, [FeFo79, BergSO, FeFo82, Ferg86, Ferg88] and 
[LLJL82, ScEu94, Cohe93, ScRo95]. Inspired by the polynomial result, but 
not the details, the first author of this paper formulated what he thought 
was a new algorithm [BaFe92, ArFe92] and gave a polynomial proof. This 
proof was independent of that of [HJLS89], a different analysis, but flawed by 
giving a slightly higher degree polynomial in the dimension than the polyno- 
mial proof given in [HJLS89]. This algorithm in [BaFe92, ArFe92] was called 
PSLQ and had the advantage of the adjustable parameter 7 or r. Applica- 
tions and implementation of this earlier version of PSLQ(r) were described 
in [BaBG94, Bail95, BaBP95]. These implementations showed that the pa- 
rameters were a helpful feature of the algorithm. The bound on iterations for 
HJLS proven in [HJLS89] was 0(n 3 +n 2 log 2 M x ); this is consonant with the 
bound proven in this paper for PSLQ(>/2). The subsequent paper [RoSc95] 
included parameters as well as addressing a certain issue of stability. 

Examples can be generated from the Mathematica implementations 
described in Section 10. Specifically, in three dimensions, consider the triple 
x = (11,27.31). We list the sequence of A -1 matrices for each algorithm. A 
relation if found will be constructed as a column of one of these A 1 matrices. 

For PSLQ( 1.1547) the successive iterations k - 0, 1, 2, 3,4, yield the five 


l-l 


matrices 



Note that PSLQ has constructed two relations appearing as the first and 
second columns of the last matrix, iteration k — 4. 

For HJLS the successive iterations fc = 0, 1, 2, 3, 4, 5, 6 yield the seven A 

matrices 
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Note that only one relation is found; it appears in the last column of the last 
matrix, iteration k = 6 . The report [ScEu 94 ] claimed that HJLS is a special 
case of PSLQ(r) for 7 = \/2 or equivalently r = yj 4 / 3 . The example just 
given shows that this claim is not true. This fact is also underscored by the 
results reported in [BaFe92], which show that HJLS, as stated in [HJLS89], 
often requires a level of numeric precision far higher than that of the input 
data, whereas PSLQ(r ) typically only requires 10 digits or so more than the 
input data. Indeed, it is clear from the results in [BaFe92] that without 
some suitable modification to the HJLS algorithm, such as that proposed in 
[ScRo95], it is not usable for many problems of interest. 

The various algorithms in the literature stand independently of any proofs. 
Though the proofs were exponential, the algorithms stated in [FeFo 79 ], and in 
[f' e f'o82], and again in [Ferg 86 ] were parametric. The parameter b in [FeFo 79 , 
FeFo82] satisfies 1 < b < 2 whereas in [Ferg 86 ] the parameter 7 is empha- 
sized. The algorithm in [Berg80, Sect, x] seems closest to PSLQ(i/ 473) 
with the r parameter set by 7 = \/2. This parameter choice appears in 
[Berg80, Sect, x] without the [LLj L82] setting and reappears in [HJLS89] as 
the small integer relation algorithm”, which we call HJLS, rewritten in the 
[LLjL82] language and accompanied by a polynomial time proof. 

Bergman discussed the complex case of finding gaussian integer relations 
for complex vectors in [Berg80, Sect. 5: Variants]. Bergman also gave an al- 
gorithm for the simultaneous real vector case in [Berg80, Sect. 7 ], Following 
Bergman, the paper defining HJLS for simultaneous real vectors, [HJLS89, 
cf., Sect. 5], implicitly includes the complex and quaternion vector case 
as well. As an alternate approach, inspired by [Shim 94 ], in this paper we 
have extended the base field of PSLQ(r) to these division rings and intro- 
duced unitary matrices into the algorithm directly. The proof given here 
of polynomial number of iterations covers the real and complex cases, but 
fails for quaternions. However, the quaternion version of PSLQ(r) performs 
reasonably well experimentally in finding hamiltonian integer relations for 
quaternion vectors. This was explained in Section 8 . 

10. Working Mathematical code for the 
algorithms PSLQ, PSOS, HJLS, BRUN. 

Attached to this paper as an appendix is a list of working Mathematica™ 
piocedures for PSLQ and a few other lattice relation algorithms. Each algo- 
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rithm algo is given by an initialization procedure initalgo[ ] followed by 
one iteration procedure algo [ ] . The input is some tuple such as 

1/ max \h JtJ {H)\,k) 

1 < j < n — 1 

where k is the iteration number, where h ]t} {H), 1 < j < n - 1, are the 
diagonal elements for a matrix A € GL(n,Z ) with xH = 0 for H € M(n x 
n - l, R)< H a rank n - 1 and n x n - 1 lower trapezoidal matrix, xH = 0. 

The output is the sextuple 

(xB~\A~ i B^\BA,BH, 1/ max )|, fc + 1) 

for the matrix A € GL(n, Z). PSLQ and PSOS are written in real GL{n , Z) 
and complex GL{n, Z -f iZ) form. Variations from this format will be clear 
from the actual Mathematical^ notation itself. 

Each algorithm is characterized by the sequence of matrices developed at 
each iteration, from the group GL(n, Z) or the group GL(n , Z + tZ), 


A u A 2 ,.-- ,A*,... 

(“partial quotients”) and their accumulated products ( convergents ) 

Bi, B 2 , . . . ,Bk-- ■ 

The general scenario is that B k H x will converge to zero if the coordinates 
of x are integrally linearly independent, B k H x gives lower bounds on the 
size of possible relations, the rows of B k will converge to x, and if x has a 
relation at all. then one will appear as a column of B k 1 for some k. There are 
counterexamples for BRUN (relations not always found), no counterexamples 
are known for PSOS. This paper proves that PSLQ(r) fits this scenario where 
k is polynomial^ bounded by (") log r (7 n ~ l M x ) with M x the least norm of 
any relation. HJLS fits this scenario as well, cf., [HJLS89]. 

In these listings, sO [x, j ,d] computes the square root of the j- th partial 
sum of squares for x, precision d decimals, x0[x,d] computes the unit vec- 
tor with the same direction as x, h0[x,d] constructs the initial H x matrix 
out of the partial sums of squares, hLQ [H] is the maximum of the absolute 
values of the diagonal elements in the LQ decomposition of H. By Theorem 
1, the reciprocal of this number gives a lower bound on the I 2 norm of any 
relation for x. brun is Brun’s algorithm: this is a simple generalization of 
Euclid’s anthyphairesis algorithm n = 2 to n > 2 quantities. Anthyphairesis, 
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or adri^cciptais, means “continually subtracted in turn from", cf.,[ Fowl 79 ], 
This concept appears in Euclid in Book VII and Book X, [Euclid], Brun’s 
algorithm, cf., [Brunl9], [Brun57], is a natural generalization of aOpcpoa pea a 
from a pair of numbers to a list of numbers. This generalization is rediscov- 
ered by almost everyone working in this area. Brun’s algorithm can cycle and 
does not always find relations. However, according to a theorem of Forcade, 
cf., [Forc81], Brun s algorithm finds relations almost everywhere, pslq is 
the PSLQ algorithm described in this paper; tau and rho are the PSLQ pa- 
rameters defined in this paper, psos is the partial sum of squares algorithm 
defined in [Ferg 88 ]. hjls is the “small integer relation algorithm” defined in 
[HJLS89]; see also [Berg80], [Ferg87]. 

11. Open Questions 

1 ) Is there a relation finding algorithm that finds one of the shortest 
relations (there may be more than one with the same minimum height), with 
a guaranteed iteration count that is a polynomial function of the dimension? 

2 ) What are the best choices for the parameter r or 7 relative to the 
number of iterations, time, and precision requirements of PSLQ? 

3) Does PSOS have a counterexample in dimension 5 or less? The com- 
plete Mathematical definition of PSOS for real and complex numbers, with 
possible extension to quaternions, is described in Section 10 . 
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